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Abstract 

The inviscid limit of the stochastic Burgers equation, with body 
forces white noise in time, is discussed in terms of the level surfaces 
of the minimising Hamilton-Jacobi function, the classical mechanical 
caustic and the Maxwell set and their algebraic pre-images under the 
classical mechanical flow map. The problem is analysed in terms of a 
reduced (one dimensional) action function. We give an explicit expres- 
sion for an algebraic surface containing the Maxwell set and caustic 
in the polynomial case. Those parts of the caustic and Maxwell set 
which are singular are characterised. We demonstrate how the geome- 
try of the caustic, level surfaces and Maxwell set can change infinitely 
rapidly causing turbulent behaviour which is stochastic in nature, and 
we determine its intermittence in terms of the recurrent behaviour of 
two processes. 



1 Introduction 



Burgers equation has been used in studying turbulence and in modelling the 
large scale structure of the universe [H [9l [28], as well as to obtain detailed 
asymptotics for stochastic Schrodinger and heat equations [TUl [HI [221 [SDl 
[311 [32] • It has also played a part in Arnol'd's work on caustics and Maslov's 
works in semiclassical quantum mechanics [31 [H [201 El]- 

We consider the stochastic viscous Burgers equation for the velocity field 
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t) G M"', where x G M"* and t > 0, 



— +K ■^)v^ = ^Av^'-WV{x)-eWkt{x)Wu v^{x, 0) = V5o(x)+0(/i'), 

Here VF^ denotes white noise and /i^ is the coefficient of viscosity which we 
assume to be small. We are interested in the advent of discontinuities in the 
inviscid limit of the Burgers fluid velocity v^{x,t) where v'^{x,t) — >■ t) 
as /i — > 0. 

Using the Hopf-Cole transformation v'^{x, t) = -/i2vin t), the Burg- 
ers equation becomes the Stratonovich heat equation, 

^ = i^Au^+fi-^V{x)u''+^kt{x)u^'oWt, u^{x, 0) = exp Ux), 
where the convergence factor Tq is related to the initial Burgers fluid density 
Now let, 

A[X]:=l- fx\s)ds- fv{X{s))ds-e [\,{Xis)) dW^, (1) 
2 Jo Jo Jo 

and select a path X which minimises ^[^]. This requires, 

dX{s) + VV{X{s)) ds + eVh{X{s)) dWs = 0. (2) 

We then define the stochastic action A{X{0),x,t) := inf {^[X] : X(t) = x} . 
Setting, 

^(X(0), X, t) := ^o(X(0)) + A(X(0), X, t), 

and then minimising A over X{0), gives X{0) = \/So{X{0)). Moreover, it 
follows that, 

St{x) := inf {^(X(0),x,t)}, 

X{0) 

is the minimal solution of the Hamilton- Jacobi equation, 

dSt + Q|V5j|' + V{x)^ dt + ekt{x) dWt = 0, St=o{x) = So{x). (3) 

Following the work of Donsker, Freidlin et al [12], — /i^ lnu^(x, t) — > St{x) 
as yU ^ 0. This gives the inviscid limit of the minimal entropy solution of 
Burgers equation as v^{x,t) = ViSf(x) [5]. 

Definition 1.1. The stochastic wavefront at time t is defined to be the set, 

Wt = {x: St{x) = 0}. 
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For small fi and fixed t, u^{x,t) switches continuously from being expo- 
nentially large to small as x crosses the wavefront Wj. However, and 
can also switch discontinuously. 

Define the classical flow map : M'^ ^ by, 

d<i>, + VV($.) ds + eWKi^,) dW^. = 0, $0 = id, $0 = V^o- 

Since X{t) = x it follows that X{s) = $s {^T^i^))^ where the pre-image 
xo{x,t) = <l>7^(x) is not necessarily unique. 

Given some regularity and boundedness, the global inverse function the- 
orem gives a caustic time T{uj) such that for < t < T{uj), $t is a random 
diffeomorphism; before the caustic time v^{x,t) = <j>t (^^^(x)) is the inviscid 
limit of a classical solution of the Burgers equation with probability one. 

The method of characteristics suggests that discontinuities in v^{x,t) are 
associated with the non-uniqueness of the real pre-image Xo{x, t). When this 
occurs, the classical flow map $t focusses an infinitesimal volume of points 
dxQ into a zero volume dX(t). 

Definition 1.2. The caustic at time t is defined to he the set, 

Assume that x has n real pre-images, 

<^^^{x] = {xo{l){x,t),xo{2){x,t), . . . ,xo{n){x,t)} , 

where each xo{i){x,t) G W^. Then the Feynman-Kac formula and Laplace's 
method in infinite dimensions give for a non-degenerate critical point, 

u^{x,t) = Y,^O.e^p(^-^^y (4) 

where Sq{x, t) := Sq {xo{i){x, t)) + A {xo{i){x, t),x, t) , and 9i is an asymptotic 
series in /i^. An asymptotic series in /i^ can also be found for v^{x,t) [SB] . 
Note that St{x) = min{S'Q(x, t) : i = 1,2, . . . , n}. 

Definition 1.3. The Hamilton- Jacobi level surface is the set, 

HI = |x : Sq{x, t) = c for some z} . 
The zero level surface H^ includes the wavefront Wt- 
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As /i ^ 0, the dominant term in the expansion (jlj) comes from the min- 
imising xo{i){x,t) which we denote Xo{x,t). Assuming xo{x,t) is unique, we 
obtain the inviscid hmit of the Burgers fluid velocity as v^{x, t) = $f (xo(x, t)) . 

If the minimising pre-image xo{x,t) suddenly changes value between two 
pre-images xo{i){x,t) and xo{j){x,t), a jump discontinuity will also occur 
in the inviscid limit of the Burgers fluid velocity. There are two distinct 
ways in which the minimiser can change; either two pre-images coalesce and 
disappear (become complex), or the minimiser switches between two pre- 
images at the same action value. The first of these occurs as x crosses the 
caustic and when the minimiser disappears the caustic is said to be cool. The 
second occurs as x crosses the Maxwell set and again, when the minimiser is 
involved the Maxwell set is said to be cool. 

Definition 1.4. The Maxwell set is given by, 

Mt = {x: 3xo,Xo G M'^ s.t. 

X = $t(xo) = $t(xo), Xq ^ Xq and A{xo, x, t) = A{xo, x,t)} . 

Example 1.5 (The generic Cusp). Let V{x, y) = 0, kt{x, y) = and Sq{xo, yo) = 
Xgyo/2. This initial condition leads to the generic Cusp, a semicubical 
parabolic caustic shown in Figure [H The caustic Ct (long dash) is given 
by, ^ ^ 

Xt{xo) = t^Xo, yt{xo) = -txl - -. 

The zero level surface (solid line) is, 

X(tfi){xQ) = y ± 1 - f^xi^ , y(t,o)(a;o) = ^ [t^xl - 1 ± y'T^t^^^ , 

and the Maxwell set Mt (short dash) is x = for y > —1/t. 

Notation: Throughout this paper x, xo,Xf etc will denote vectors, where 
normally x = $((xo). Cartesian coordinates of these will be indicated 
using a sub/superscript where relevant; thus x = {xi,X2, ■ ■ ■ ,Xd), Xq = 
(xq, Xq, . . . , Xq) etc. The only exception will be in discussions of explicit 
examples in two and three dimensions when we will use {x,y) and (xo,?/o) 
etc to denote the vectors. 

2 Some background 

We begin by summarising some of the geometrical results established by 
Davies, Truman and Zhao (DTZ) P, [3 [H] and presenting some minor gener- 
alisations of their results [221 [25]. Following equation ([T]), let the stochastic 
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Figure 1: Cusp and Tricorn. 



action be defined, 

Aixo,Po,t) = l- [ X{s)'ds- [ \viXis))ds + ek,iXis))dW, 
Jo Jo . 

where X{s) = X{s,Xo,po) G M'^ and, 

dX{s) = -VViX{s))ds-eVh{X{s))dWs, X(0) = xq, X{0) = po, 

for s G [0, t] witli Xo,po G M.'^. We assume X{s) is JF^ measurable and unique. 

Lemma 2.1. Assume So,V G and kt G C^'^ , W, V/c* Lipschitz with 
Hessians V'^V, V'^kt and all second derivatives with respect to space variables 
of V and kf bounded. Then for pq, possibly xq dependent, 

^(xo, Po, t) = X{t) ■ ^ - X„(0), a = 1, 2, . . . , rf. 



Metliods of Kolokoltsov et al [T8J [19] guarantee that for small t the map 
Po ^ ^it,XQ,po) is onto for all Xq. Therefore, we can define, 

A{xo,x,t) = ^(a;o,Po,^)U=po(,o,,,t), 

where po = po{xo, x, t) is the random minimiser (assumed unique) of A{xo,po, t) 
when X{t, Xo,Pq) = x. The stochastic action corresponding to the initial mo- 
mentum VS'o(xo) is then A{xo,x,t) := A{xo,x,t) + Sq^Xq). 

Theorem 2.2. If^t is the stochastic flow map, then $t(xo) = x is equivalent 
to, 

d 

-^[A{xo,x,t)] = 0, a = l,2,...,d. 
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The Hamilton- Jacobi level surface is obtained by eliminating xq be- 
tween, 

dA 

A{xo,x,t)=c and ——^{xo,x,t)=0, a = 1,2, . . . ,d. 

OXq 

Alternatively, if we eliminate x to give an expression in xq, we have the pre- 
level surface ^^^Hj:. Similarly the caustic Ct (and pre-caustic <l>^^Ct) are 
obtained by eliminating xq (or x) between, 

■^^M — ^(^o,a;,t) =0 and —^(xq, x, t) = a = l,2,...,d. 

Kdx^dx^ /a,/3=l,2,...,d 

These pre-images are calculated algebraically which are not necessarily the 
topological inverse images of the surfaces Ct and under $j. 

Assume that A{xo,x,t) is in space variables with det ^^^^j 7^ 0. 

Definition 2.3. A curve x = x(7), 7 G A^(7o, S), is said to have a generalised 
cusp at'j = 7o, 7 being an intrinsic variable such as arc length, if ■^{■Jo) = 0. 

Lemma 2.4. Let $4 denote the flow map and let ^^^Tt and Tt be some sur- 
faces where if Xq E ^^^Tt then x = $t(xo) G Tt- Then $4 is a differentiable 
map from ^t^^t to Tt with Frechet derivative, 

Lemma 2.5. Let Xq{s) be any two dimensional intrinsically parameterised 
curve, and define x{s) = $t(xo(s)). Let Cq denote the zero eigenvector of 

(dx"^^ ) ^'^^ assume that ker ^ (dl^^ ) ~ ^^o)' Then, there is a generalised 
cusp on x{s) when s = a if and only if either: 

1. there is a generalised cusp on Xq{s) when s = a; or, 

2. a;o(cr) is on the pre-caustic and the tangent -^(s) at s = a is parallel 

to Cq. 

Proposition 2.6. The normal to (^^^H^ is. 
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Corollary 2.7. In two dimensions, let $^ iJf meet $^ Ct at Xq where 
n{xo) 7^ and ker ^^^^j = (cq). Then the tangent to ^^^Hf at xq is 
parallel to Cq. 

Proposition 2.8. In two dimensions, assume that n{xo) ^ where xq G 
^^^H^, so that does not have a generalised cusp at xq. Then 

can only have a generalised cusp at $t(xo) if $t(xo) G Ct- Moreover, if 
X = $t(xo) € $t n Hf^^ then will have a generalised cusp. 

Example 2.9 (The generic Cusp). Figure [2] sliows tliat a point lying on 
three level surfaces has three distinct real pre-images each on a separate pre- 
level surface. A cusp only occurs on the corresponding level surface when 
the pre-level surface intersects the pre-caustic. Thus, a level surface only has 
a cusp on the caustic, but it does not have to be cusped when it meets the 
caustic. 




Figure 2: (a) The pre-level surface (solid line) and pre-caustic (dashed), (b) 
the level surface (solid line) and caustic (dashed), both for the generic Cusp 
with c > 0. 



Theorem 2.10. Let, 

X e Cusp{H^) = {xe^t {^;'Ct n ^;'H^) ,x = $t(xo),n(xo) ^ O} . 

Then in three dimensions in the stochastic case, with probability one, T^ the 
tangent space to the level surface at x is at most one dimensional. 
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3 A one dimensional analysis 



In this section we outline a one dimensional analysis first described by Reynolds, 
Truman and Williams (RTW) [H]. 

Definition 3.1. The d- dimensional flow map is globally reducible if for 
any X = {xi,X2, ■ ■ ■ ,Xd) and xo = (xq, Xq, . . . , Xq) where x = ^^(xo), it is 
possible to write each coordinate Xq as a function of the lower coordinates. 
That is, 

X = $t(xo) Xq = Xq{x,xI,xI, . . .XQ~^,t) for a = d,d- 1, . . . ,2. (5) 

Therefore, using Theorem 12.2^ the flow map is globally reducible if we 
can find a chain of functions Xq, Xq~^, . . . , Xq such that, 

— rpd,/ 1 2 ^d-l j.\ 
•^0 — •^OV'^' •^0' •^0' • • • -^0 ' 

^d-l _ ^(1-1/ 1 2 d-2 ,\ 

•^0 — ■^0 l-^; -^0' •^0' • • • "^0 ' 

^^(Xq, Xq, Xq(x, Xq, Xq, t) Xq (...), X, t) = 0, 

where Xq(. . .) is the expression only involving Xq and Xq gained by substituting 
each of the functions Xq, . . . , Xq~^ repeatedly into Xq(x, Xq, Xq, . . . , Xq~ \ t). 
This requires that no roots are repeated to ensure that none of the second 
derivatives of A vanish. We assume also that there is a favoured ordering of 
coordinates and a corresponding decomposition of which allows the non- 
uniqueness to be reduced to the level of the Xq coordinate. This assumption 
appears to be quite restrictive. However, local reducibility at x follows from 
the implicit function theorem and some mild assumptions on the derivatives 
of A. 

Definition 3.2. // is globally reducible then the reduced action function 
is the univariate function gained from evaluating the action with equations 

f{x,t){xQ) := /(Xq, X, t) = A{Xq, Xq(x, Xq, t), Xq(. ..),... ,x,t). 



^ 77-7(3:0, a;, t) = 0, 



dx[^ 
OA 



I (a^o, ^0' • • • 5 a;Q(. . .), X, t) — 0, 
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Lemma 3.3. // $t is globally reducible, modulo the above assumptions, 



clet I — — —{XQ,x,t) 



\{dxoy 



n 



^0 — (^0''^of'^''^0'^^'*"'''^of '■'^■^ 



d 



a ] -^{^Qi ■ ■ ■ i^Oi^O^ {■■■)■,■■■ iXq{- ■ ■)-iX,t) 



where the first term is f'^^f^ixl) and the last d — 1 terms are non zero. 

Theorem 3.4. Let the classical mechanical flow map $t be globally reducible. 
Then: 

1. f[xf^{x\) = and the equations 1^ -v^ x = $f(xo), 

^- f(x,t)(^o) = f(x,t)i^o) = ^'^^ ^^'^ equations ^ 

X = $t(xo) is such that the number of real solutions Xq changes. 



4 Analysis of the caustic 

We begin by parameterising the caustic = det (D$((xo)) from Definintion 
II. 2t this equation only involves Xq and t, and is therefore the pre-caustic. We 
use this to parameterise the pre-caustic as, 

Xq = Al, Xq = A2, Xq ^ = Xd-l aud Xq = Xq ( A 1 , A2 , • • • , Arf. 1 ) . 

The parameters are restricted to be real so that only real pre-images are 
considered. 

Definition 4.1. For any A = (Ai, A2, . . . , \d~i) G M'^^^ the pre-parameterisation 
of the caustic is given by Xt{\) := (A,Xq(A)) . 

The pre-parameterisation will be intrinsic if ker(Z}$t) is one dimensional. 

Corollary 4.2. Let Xt{X) denote the pre-parameterisation of the caustic 
where A = (Ai, A2, . . . , A,„i) G R^'K Then 4,(,),i)(Ai) = /(';,(a),*)(Ai) = 0. 

Proposition 4.3. Let Xt{X) denote the pre-parameterisation of the caustic 
where A = (Ai, A2, . . . , Ad_i) G IR'^"^ Assume f{xt{x),t){xl) G C^~^^ then, in 
d-dimensions, if the tangent to the caustic is at most {d — p+V}- dimensional 
at Xt{\), 

4,(A),t)(^l) = /(",(A),t)^^l) = • • • = /(t\A),t)^^l) = 0- 
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Proof. Follows by repeatedly differentiating f'(xt{\)t)^^^) ~ 0' "which holds if 
the tangent space at Xt(A) is {d — 2)-dimensional [22]. □ 

From Corollary 14.21 there is a critical point of inflexion on f{x,t){xl) at 
Xg = Ai when x = Xt{X). Consider an example where for x on one side of 
the caustic there are four real critical points on f\x,t){xo) = 0. Let them be 
enumerated XQ{i){x, t) for i = 1 to 4 and denote the minimising critical point 
xl{x,t). Figure [H] illustrates how the minimiser jumps from (a) to (6) as x 
crosses the caustic. This will cause u'^ and v'^ to jump for small /i and the 
caustic at such a point is described as being cool. 



Before Caustic 



On Cool Caustic Beyond Caustic 




!,(2T=x'i(x,i) 

Minimiser at 
xli2){x,t) = xlix,t). 





Two 's coalescing 
form point of inflexion. 



Minimiser jumps. 



Figure 3: The graph of f{x,t){xl) as x crosses the caustic. 



Definition 4.4. Let Xt{\) be the pre-parameterisation of the caustic. Then 
Xt{X) IS on the cool part of the caustic if f{xt{x),t){>^i) < f{xt{\),t){xl{i){xt{\),t)) 
for all i = 1,2, ... ,n where xl{i){x, t) denotes an enumeration of all the real 
roots for Xq to f[xt)i.^})) = 0- // caustic is not cool it is hot. 

Definition 4.5. The pre-normalised reduced action function evaluated on 
the caustic is given by J^x{xl) := f(xtix),t){xl) - /(xt(A),t)(Ai). 

Assume that J-'x( real analytic function in a neighbourhood of 

Ai e R. Then, 

J'x{xl) = {xl-\,fF{xl), (6) 

where F is real analytic. When the inflexion at Xq = Ai is the minimis- 
ing critical point of T\, the caustic will be cool. Therefore, on a hot/cool 
boundary this inflexion is about to become or cease being the minimiser. 

Proposition 4.6. A necessary condition for Xt{X) G Ct to be on a hot/cool 
boundary is that either F{xl) or G{xl) has a repeated root at Xq = r where, 

G{xl) = ?,F{xl) + {xl-\{)F'{xl). 
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Proof. The minimiser could change when either F has a repeated root which 
is the minimiser, or there is a second inflexion at a lower minimising value 
[23]. □ 

The condition is not sufficient as it includes cases where the minimiser is 
not about to change (see Figure Hj) . 



Increasing A 



Possible 
hot/cool 
boundary 



Caustic changes hot to cool 




No change in caustic 




Figure 4: Graphs of J^\{xl) as A varies. 



Example 4.7 (The polynomial swallowtail). Let V{x,y) = 0, kt{x,y) = x, 
and 5'o(a;o,2/o) = + x^yo. This gives global reducibility and kt{x,y) = x 
means that the effect of the noise is to translate e = picture through 
^— e Ws ds, 0^ . A simple calculation gives, 

F{xo) = 12A2 - 3Xt + QXxo - txo + 2x1, 
G{xo) = ISA^ - 4At + lOAxo - 2txo + Sx^. 




Hot 



Figure 5: Hot and cool parts of the polynomial swallowtail caustic for t = 1. 
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Example 4.8 (The three dimensional polynomial swallowtail). Let V{x, y) = 
0, kt{x, y) = 0, and 50(2:0, yo, Zq) = Xq + x^yo + x^zq. The functions F and G 
can be easily found, and an exact expression for the boundary extracted [22] 
this is shown in Figure El 




Boundary on the caustic. Hot and cool parts. 



Figure 6: The hot (plain) and cool (mesh) parts of the 3D polynomial swal- 
lowtail caustic at time t = 1. 



5 Swallowtail perestroikas 

The geometry of a caustic or wavefront can suddenly change with singularities 
appearing and disappearing [2]. We consider the formation or collapse of 
a swallowtail using some earlier works of Cayley and Klein. This section 
provides a summary of results from [23] where all proofs can be found. 

We begin by recalling the classification of double points of a two dimen- 
sional algebraic curve as acnodes, crunodes and cusps (Figure [7j). 




Acnode. Crunode. Cusp. 



Figure 7: The classification of double points. 

In Cayley's work on plane algebraic curves, he describes the possible triple 
points of a curve [27] by considering the collapse of systems of double points 
which would lead to the existence of three tangents at a point. The four 
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possibilities are shown in Figure [HI The systems will collapse to form a triple 
point with respectively, three real distinct tangents, three real tangents with 
two coincident, three real tangents all of which are coincident, or one real 
tangent and two complex tangents. It is the interchange between the last 
two cases which will lead to the formation of a swallowtail on a curve fL5\ . 
This interchange was investigated by Felix Klein [T7]. 






Figure 8: Cayley's triple points. 

In Section 3, we restricted the pre-parameter to be real to only consider 
points with real pre-images. This does not allow there to be any isolated 
double points. We now allow the parameter to vary throughout the complex 
plane and consider when this maps to real points. We begin by working with 
a general curve of the form x(A) = (xi(A), X2(A)) where each Xa{X) is real 
analytic in A G C. If lm{x{a + i?])} = 0, it follows that x{a + it]) = x{a — it]), 
so this is a "complex double point" of the curve x(A). 

Lemma 5.1. If x{X) = (xi(A), X2(A)) is a real analytic parameterisation of 
a curve and A is an intrinsic parameter, then there is a generalised cusp at 
X = Xq if and only if the curves, 

= -Im{xa(a + irj)} a = 1,2, 
V 



intersect at (Ao,0) in the {a,ri) plane. 

Now consider a family of parameterised curves Xt{X) = (x|(A), x^(A)). As 
t varies the geometry of the curve can change with swallowtails forming and 
disappearing. 

Proposition 5.2. // a swallowtail on the curve Xt{X) collapses to a point 
where A = A when t = i then, 

"^"^^-(A) = ^(A) = 0. 



dA' ' dA2 
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Proposition 5.3. Assume that there exists a neighbourhood of X E M. such 
that -^(A) 7^ fort E (t — 5,i) where S > 0. If a complex double point joins 
the curve Xt{X) at X = X when t = t then, 

These provide a necessary condition for tlie formation or destruction of a 
swallowtail, and for complex double points to join or leave the main curve. 

Definition 5.4. A family of parameterised curves Xt{X), (where X is some 
intrinsic parameter) for which, 

is said to have a point of swallowtail perestroika when X = X and t = t. 

As with generalised cusps, we have not ruled out further degeneracy at 
these points. Moreover, as Cayley highlighted, these points are not cusped 
and are barely distinguishable from an ordinary point of the curve 



5.1 The complex caustic in two dimensions 

The complex caustic is the complete caustic found by allowing the parameter 
A in the pre-parameterisation Xj(A) G to vary over the complex plane. By 
considering the complex caustic, we are determining solutions a = at and 
V = Vt to, 

+ i^) = /{x,t)(a + i^) = 0, 

where x e M^. We are interested in these points if they join the main caustic 
at some finite critical time t. That is, there exists a finite value t > such 
that rjt as t ^ i. If this holds then a swallowtail can develop at the 
critical time t. 

Theorem 5.5. For a two dimensional caustic, assume thatxt{X) is a real an- 
alytic function. If at a time t a swallowtail perestroika occurs on the caustic, 
then X = X(-(A) is a real solution for x to, 

where X = a^. 
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Theorem 5.6. For a two dimensional caustic, assume that Xt{X) is a real 
analytic function. If at a time t there is a real solution for x to, 

(A) = (A) = /(:,-)(A) = fl2,^iX) = 0, 

and the vectors Va;/(^£)(A) and Va;/('^£)(A) are linearly independent, then x 
is a point of swallowtail perestroika on the caustic. 

Example 5.7. Let V{x, y) = 0, kt{x, y) = and 5'o(xo, yo) = Xq + a^o^/o- The 
caustic has no cusps for times t < i and two cusps for times t > i where 
i = 4^2 X 33^/^ X 7(-^/^) = 2.5854 . . . 

At the critical time t the caustic has a point of swallowtail perestroika as 
shown in Figures [9] and [lOl The conjugate pairs of intersections of the curves 
in Figure M are the complex double points. There are five before the critical 
time and four afterwards. The remaining complex double points do not join 
the main caustic and so do not influence its behaviour for real times. 




Figure 9: lm{xt{a + irj)} = (sohd) and lm{yt{a + irj)} = (dashed) in 
(a, 77) plane. 




Figure 10: Caustic plotted at corresponding times. 



5.2 Level surfaces 

Unsurprisingly, these phenomena are not restricted to caustics. There is an 
interplay between the level surfaces and the caustics, characterised by their 
pre-images. 
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Proposition 5.8. Assume that in two dimensions at xq G $7 fl $7 

the normal to the pre-level surface ri(xo) 7^ and the normal to the pre- 
caustic n{xo) ^ so that the pre-caustic is not cusped at xq- Then n{xo) is 
parallel to n{xo) if and only if there is a generalised cusp on the caustic. 

Corollary 5.9. Assume that in two dimensions at Xq G ^^^H^ fl ^i^^Ct the 
normal to the pre-level surface n{xo) 7^ 0. Then at $t(xo) there is a point 
of swallowtail perestroika on the level surface if and only if there is a 
generalised cusp on the caustic Ct at $((xo). 

Example 5.10. Let V{x,y) = 0, kt{x,y) = 0, and So{xo,yo) = x^ + x^yo. 
Consider the behaviour of the level surfaces through a point inside the caustic 
swallowtail at a fixed time as the point is moved through a cusp on the 
caustic. This is illustrated in Figure [TH Part (a) shows all five of the level 
surfaces through the point demonstrating how three swallowtail level surfaces 
collapse together at the cusp to form a single level surface with a point of 
swallowtail perestroika. Parts (b) and (c) show how one of these swallowtails 
collapses on its own and how its pre-image behaves. 




Figure 11: (a) All level surfaces (solid line) through a point as it crosses 
the caustic (dashed line) at a cusp, (b) one of these level surfaces with its 
complex double point, and (c) its real pre-image. 



6 Maxwell sets 

A jump will occur in the inviscid limit of the Burgers velocity field if we 
cross a point at which there are two different global minimisers xo{i){x,t) 
and xo{j){x,t) returning the same value of the action. 
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In terms of the reduced action function, the Maxwell set corresponds to 
values of x for which f(x,t){xl) has two critical points at the same height. If 
this occurs at the minimising value then the Burgers fluid velocity will jump 
as shown in Figure [121 

Before Maxwell set On Cool Maxwell set Beyond McLxwell set 




Minimiser at xj. Two xq 's at same level. Minimis er jumps. 



Figure 12: The graph of the reduced action function as x crosses the Maxwell 
set. 

6.1 The Maxwell-Klein set 

We begin with the two dimensionals polynomial case by considering the clas- 
sification of double points of a curve (Figure [7]). 

Lemma 6.1. A point x is in the Maxwell set if and only if there is a 
Hamilton- Jacobi level surface with a point of self-intersection (crunode) at 

X. 

Proof. Follows from Definition II. 4[ □ 

Definition 6.2. The Maxwell-Klein set Bf is the set of points which are 
non-cusp double points of some Hamilton- Jacobi level surface curve. 

It follows from this definition that a point is in the Maxwell-Klein set 
if it is either a complex double point (acnode) or point of self-intersection 
(crunode) of some Hamilton- Jacobi level surface. Using the geometric results 
of DTZ outlined in Section 2, it is easy to calculate this set in the polynomial 
case as the cusps of the level surfaces sweep out the caustic. 

Theorem 6.3. Let Dt be the set of double points of the Hamilton- Jacobi 
level surfaces, Ct the caustic set, and Bt the Maxwell-Klein set. Then, from 
Cayley and Klein's classification of double points as crunodes, acnodes, and 
cusps, by definition Df = Ct U Bt and the corresponding defining algebraic 
equations factorise Dt = ■ B^, where m,n are positive integers. 
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Proof. Follows from Proposition 12.81 and Lemma 16. 1[ □ 
Theorem 6.4. Let p(t,c){x) be the resultant, 

P(t,c){x) = R - c, /(^,t)(-)) , 

where x = [xi, X2)- Then x E Dt if and only if for some c, 

P(t,c)(x)-^(x)-^(x)-0. 

Further, 

Dt{x) = gcd (pj(x),p?(x)) , 

where gcd(-,-) denotes the greatest common divisor and p\ and pi are the 
resultants, 

p]{x) = R [p,,,{xl ^(x)) and p]{x) = R (^(x), ^-^{x, 

Proof. Recall that the equation of the level surface of Hamilton- Jacobi func- 
tions is merely the result of eliminating x\ between the equations, 

f{x,t){.xl) = c and f[^^t)i.xl) = 0. 

We form the resultant p(^t,c) (x) using Sylvester's formula. The double points 
of the level surface must satisfy for some c G M, 

P,,e)(x) = 0, ^(x) = and (x) = 0. 

Sylvester's formula proves all three equations are polynomial in c. To proceed 
we eliminate c between pairs of these equations using resultants giving, 



^(p(..)(-),^(-))=pK-) and r(^ 



„2, 



Let Dt = gcd(pj , pf) be the greatest common divisor of the algebraic p] and 
pI- Then Dt{x) = is the equation of double points. □ 

We now extend this to (i-dimensions, where the Maxwell-Klein set corre- 
sponds to points which satisfy the Maxwell set condition but have both real 
pre-images (Maxwell) or complex pre- images (Klein). 
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Theorem 6.5. Let the reduced action function f{x,t){xo) be a polynomial 
in all space variables. Then the set of all possible discontinuities for a d- 
dimensional Burgers fluid velocity field in the inviscid limit is the double 
discriminant, 

Dit):=D,{Dx (/(,,i)(A)-c)} =0, 
where Dx{p{x)) is the discriminant of the polynomial p with respect to x. 

Proof. By considering the Sylvester matrix of the first discriminant, 

m 

Dx (/(.,t)(A) -c)=K\[ - c) , 

i=l 

where XQ{i) (x, t) is an enumeration of the real and complex roots A of /(^. (A) = 
and K is some constant. Then the second discriminant is simply, 

{D, (/(,,,)(A) - c)) = K'-^-' n ifi.d^i^M) - f^,,){xl{j){x,t))Y . 

□ 

Theorem 6.6. The double discriminant D{t) factorises as, 

D{t) = bl^-'-{Ctf-iBt)\ 

where Bf = is the equation of the Maxwell-Klein set and Ct = is the 
equation of the caustic. The expressions Bt and Ct are both algebraic in x 
and t. 

Proof See [23]. □ 

Example 6.7 (The polynomial swallowtail). Let V{x,y) = 0, kt{x,y) = 
and, S'o(xo, yo) = ^1 + XqUo- The Maxwell-Klein set can be found by factori- 
sation giving, 

= -675 + 52t^-t* + 3120t^x- 224t^x + 4t^^x- 38400tV + 1408tV 
+ 128000tx^ - 5400ty + 312t^y - At^y + l2Amt^xy - AA^t^xy 
-reSOOt^x^y - 16200t2?/2 + Q2Afy'^ - At^^y"^ + UAm^xy"^ 
-21600t^y^ + 416tV - lOSOOtV- 

Outside of the swallowtail on the caustic there are two real and two complex 
pre-images whereas inside there are four real and no complex pre-images. 
Therefore, any part of the Maxwell-Klein set outside of the caustic swallowtail 
must correspond to Klein double points and any part inside must correspond 
to the Maxwell set. This is shown in Figure [T3l 
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■Acnodes of 

= Klein set. 



Figure 13: The caustic and Maxwell-Klein set. 
6.2 The pre-Maxwell set 

If the Maxwell set is defined as in Definition II. 4[ then the pre-Maxwell set is 
the set of all the pre-images Xq and Xq which give rise to the Maxwell set. 

Definition 6.8. The pre-Maxwell set ^^^Mf is the set of all points Xq G M'^ 
where there exists x,xq G M.'^ such that x = $i(xo) and x = $i(xo) with 
Xq ^ Xq and, 

A{xo, X, t) = A{xo, X, t). 

With the caustic and level surfaces, each regular point was linked by 
to a single point on the relevant pre-surafce. However, every point on the 
Maxwell set is linked by to at least two points on the pre-Maxwell set. 

Theorem 6.9. The pre-Maxwell set is given by the discriminant D^i {G{xl)) = 
where, 



G{xl) 



[x^-xiy 



Proof. From the Definition 16.81 and Theorem 13.41 it follows that the pre- 
Maxwell set is found by eliminating x and x^ between, 

f{x,t){^l) = f{x,t){^l) fl^,t){xl) = /(^,j)(Xo) = 0- 

This surface would include the pre-caustic where Xg = Xg and so this repeated 
root must be eliminated. □ 
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We can use this to pre-parameterise the Maxwell set as has been done 
with the caustic and level surfaces. By restricting the parameter to be real, 
we only get the Maxwell set as the Klein points have complex pre-images. 

We now summarise the results of 1251. 



Lemma 6.10. Assume that a point x on the Maxwell set corresponds to 
exactly two pre-images on the pre-Maxwell set, xq and xq. Then the normal 
to the pre-Maxwell set at Xq is to within a scalar multiplier given by, 

X(t, xo, VS'o(xo)) - X{t, xo, VS'o(xo)) 



Corollary 6.11. In two dimensions let the pre-Maxwell set meet the pre- 
caustic at a point Xq where n 7^ and 

( d^A \ 
ker ( y^(xo, $i(xo),t) 1 = (cq), 

where cq is the zero eigenvector. Then the tangent plane to the pre-Maxwell 
set at Xq, Tj-o is spanned by Cq. 

Proposition 6.12. Assume that in two dimensions at xq G ^^^Mt the nor- 
mal n(xo) ^ so that the pre-Maxwell set does not have a generalised cusp 
at Xq. Then the Maxwell set can only have a cusp at $t(xo) «/$f(xo) G Cf. 
Moreover, if 

X = $i(xo) G $i {<l>,-'a n $r'Mi} , 
the Maxwell set will have a generalised cusp at x. 

Corollary 6.13. In two dimensions, if the pre-Maxwell set intersects the 
pre-caustic at a point xq, so that there is a cusp on the Maxwell set at the 
corresponding point where it intersects the caustic, then the pre-Maxwell set 
touches the pre-level surface ^1~^H^ at the point xq. Moreover, if the cusp 
on the Maxwell set intersects the caustic at a regular point of the caustic, 
then there will be a cusp on the pre-Maxwell set which also meets the same 
pre-level surface ^^^H^ at another point xq. 

Corollary 6.14. When the pre-Maxwell set touches the pre-caustic and pre- 
level surface, the Maxwell set intersects a cusp on the caustic. 

Example 6.15 (The polynomial swallowtail). Let V{x,y) = 0, kt{x,y) = 
and, 5*0 (xo, yo) = xl + x^yo- 
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Figure 14: The caustic (dashed) and Maxwell set (solid line). 



From Proposition 16.121 the cusps on the Maxwell set correspond to the 
intersections of the pre-curves (points 3 and 6 on Figure [T4l) . But from 
Corollary 16.131 the cusps on the Maxwell set also correspond to the cusps 
on the pre-Maxwell set (points 2 and 5 on Figure [T3] and also Figure [T5l) . 
The Maxwell set terminates when it reaches the cusps on the caustic. These 
points satisfy the condition for a generalised cusp but, instead of appearing 
cusped, the curve stops and the parameterisation begins again in the sense 
that it maps back exactly on itself. At such points the pre-surfaces all touch 
(Figure fT5l) . 




Cusp on Maxwell set Cusp on caustic 



Figure 15: The caustic (long dash) and Maxwell set (solid line) with the level 
surfaces (short dash) through special points. 

These two different forms of cusps correspond to very different geometric 
behaviours of the level surfaces. Where the Maxwell set stops or cusps corre- 
sponds to the disappearance of a point of self-intersection on a level surface. 
There are two distinct ways in which this can happen. Firstly, the level sur- 
face will have a point of swallowtail perestroika when it meets a cusp on the 
caustic. At such a point only one point of self-intersection will disappear. 
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and so there will be only one path of the Maxwell set which will terminate 
at that point. However, when we approach the caustic at a regular point, 
the level surface must have a cusp but not a swallowtail perestoika. This 
corresponds to the collapse of the second system of double points in Figure 
[HI Thus, two different points of self intersection coalesce and so two paths of 
the Maxwell set must approach the point and produce the cusp (see Figure 




Approaching the caustic Approaching a cusp on the caustic 



Figure 16: The caustic (long dash) Maxwell set (solid line) and level surface 
(short dash). 



7 Some applications to turbulence in two di- 
mensions 

7.1 Real turbulence and the ( process 

Definition 7.1. The turbulent times t are times when the pre-level surface 
of the minimising Hamilton- J acohi function touches the pre-caustic. Such 
times t are zeros of a stochastic process C'^(-)- ^-S- C'^(^) = 0. 

These turbulent times are times at which the number of cusps on the 
corresponding level surface will change. We begin with some minor general- 
isations of results in RTW [3l] and also [231 126] . 

Proposition 7.2. Assume is globally reducible and that Xt{X) is the pre- 
parameterisation of a two dimensional caustic. Then the turbulence process 
at A is given by, 

Cit) = /(xt{Ao),t)(Ao) - c, 

where f{x,t){,x\) is the reduced action evaluated at points x = Xt(Ao) where 
Xt{\o) = $t(Ao,Xo(Ao)) e Ct, A = Ao satisfying, 

^.(A)-^(A) = 0, 
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where Xt{X) = i't(A, Xq q{X)) and Xt{Xo) G , the cool part of the caustic. 
Hence, there are three kinds of real stochastic turbulence:- 

1. Cusped, where there is a cusp on the caustic, 

2. Zero speed, where the Burgers fluid velocity is zero, 

3. Orthogonal, where the Burgers fluid velocity is orthogonal to the caus- 
tic. 

Proof. The number of cusps on the relevant pre-level surface is, 

n,(t) = #{AGM:/(.,(A),t)(A) = c}, 

where the roots A = Aq correspond to points in the cool part of the caustic. 
The pre-surfaces touch when ndt) changes, which occurs when, 

^/(x.(A),)(A)=0. □ 

For stochastic turbulence to be intermittent we require that the process 
C^{t) is recurrent. 

Proposition 7.3. Let V{x,y) = 0, kt{x,y) = x and 

Soixo,yo) = f{xo) + g{xo)yo, 

where f,g, f and g' are zero at xq = a but g"{a) ^ 0. Then, for orthogonal 
turbulence at a, 

ft 2 rt 

C(t) = -aeWt + e^Wt / W,ds-- ds - c. 

Jo ^ Jo 

We note the following result of RTW [Mj . 

Lemma 7.4. Let Wt be a i?M(M) process starting at 0, c any real constant 
and 

rt 2 rt 

Yt = -aeWt + e^Wt / Wsds - - ds - c. 

Jo ^ Jo 

Then, with probability one, there exists a sequence of times tn oo such that 

Yt^ = for every n. 

We also note that this can be extended to a (i-dimensional setting where 
for a (i- dimensional Wiener process W{t) the zeta process can be found ex- 
plicitly [22]. 
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Theorem 7.5. In d- dimensions, the zeta process is given by, 

Ct = f^.ow,t)i>^i)-ex',{X)-W{t) + e'W{t)- J^i^^ 

where f^^f^ixl) denotes the deterministic reduced action function, x°(A) de- 
notes the pre-parameterisation of the deterministic caustic and A must satisfy 
the equation, 

Va (/(°.o(,),)(Ai)-6a;?(A)-iy(t))=0. 

When A is deterministic, the recurrence of this process can be shown 
using the same argument as for the two dimensional case (further resuhs on 
recurrence can be found in [21]). Here we recapitulate our belief that cusped 
turbulence will be the most important. As we have shown, when the cusp 
on the caustic passes through a level surface, it forces a swallowtail to form 
on the level surface. The points of self intersection of this swallowtail form 
the Maxwell set. 

7.2 Complex turbulence and the resultant 77 process 

We now consider a completely different approach to turbulence. Let (A, Xq ^(A)) 
denote the parameterisation of the pre-caustic at time t. When, 

Zt = Im {$f(a + ir],xl (j{a + ir]))} , 

is random, the values of r]{t) for which Zt = will form a stochastic process. 
The zeros of this new process will correspond to points at which the real pre- 
caustic touches the complex pre-caustic. The points at which these surfaces 
touch correspond to swallowtail perestroikas on the caustic. When such a 
perestroika occurs there is a solution of the equations, 

fUW = f'UW = /(L't)(A) = /S)(A) = 0. 

Assuming that /(x,t)(a^o) polynomial in Xq we can use the resultant to state 
explicit conditions for which this holds [23] . 

Lemma 7.6. Let g and h be polynomials of degrees m and n respectively 
with no common roots or zeros. Let f = gh be the product polynomial. Then 
the resultant, 

R(f,n = (-1)-" (^ ^<4")';!°)(o) )" 

where N = m + n and R{g, h) ^ 0. 
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Since f[xt{x)t)(-^o) ^ polynomial in xq with real coefficients, its zeros 
are real or occur in complex conjugate pairs. Of the real roots, = A is 
repeated. So, 

/(xt(A),i)(4) = i^o - >^yQ{\t){xl)H(^x,t){xl), 
where Q is the product of quadratic factors, 

Qix,t){4) = fl{i4-^'t)' + iriir}, 

1=1 

and H(^x,t){xo) the product of real factors corresponding to real zeros. This 
gives, 

/(i;(A),)K)L-, = 2 n {(A - aD' + ivlf] i/(A,)(A). 
" 1=1 

We now assume that the real roots of H are distinct as are the complex 
roots of Q. Denoting /(1'^(a) t)(2^o) ^ by ft'{\) etc, a simple calculation 
gives 



^A(/r(A),/i'^(A))| = 

k=l j=/=k 



x\Rx{H,H')\\RxiQ,H: 



|2 



Kt being a positive constant. Thus, the condition for a swallowtail pere- 
stroika to occur is that 



0, 



where we call Prj{t) the resultant eta process. 

When the zeros of p,j(t) form a perfect set, swallowtails will spontaneously 
appear and disappear on the caustic infinitely rapidly. As they do so, the 
geometry of the cool part of the caustic will rapidly change as the A shaped 
sections typical of a swallowtail caustic appear and disappear. Moreover, 
Maxwell sets will be created and destroyed with each swallowtail that forms 
and vanishes adding to the turbulent nature of the solution in these regions. 
We call this 'complex turbulence' occurring at the turbulent times which are 
the zeros of the resultant eta process. 
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Complex turbulence can be seen as a special case of real turbulence which 
occurs at specific generalised cusps of the caustic. Recall that when a swal- 
lowtail perestroika occurs on a curve, it also satisfies the conditions for having 
a generalised cusp. Thus, the zeros of the resultant eta process must coin- 
cide with some of the zeros of the zeta process for certain forms of cusped 
turbulence. At points where the complex and real pre-caustic touch, the 
real pre-caustic and pre-level surface touch in a particular manner (a double 
touch) since at such a point two swallowtail perestroikas on the level surface 
have coalesced. 

Thus, our separation of complex turbulence from real turbulence can be 
seen as an alternative form of categorisation to that outlined in Section 7.1 
which could be extended to include other perestroikas. 
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